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We study the semileptonic decays of the lowest lying double heavy baryons using 
the relativistic three-quark model. We do not employ a heavy quark mass expansion 
but keep the masses of the heavy quarks and baryons finite. We calculate all relevant 
form factors and decay rates. 
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The semileptonic decays of heavy mesons and baryons are ideally suited to ex- 
tract the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements. The heaviest 
flavored bottom-charm 5 c -meson was observed by the CDF Collaboration [1] 
in the analysis of the decay mode B c — > J/iplv. The discovery of the S c -meson 
raises hopes that double heavy flavored baryons will also be discovered in the 
near future. The theoretical treatment of the systems with two heavy quarks is 
complicated by the fact that one cannot make use of an expansion in terms of 
the inverse heavy quark masses. Previously, nonrelativistic potential models, 
diquark approximation, QCD sum rules and nonrelativistic QCD have been 
used to describe the spectroscopy of double heavy baryons and to estimate 
the inclusive and some exclusive decay modes of such systems (for review, see 
[2]- [4] and references therein). 

In [5] we have studied the semileptonic decays of the double heavy £> c -meson 
within a relativistic constituent quark model. The relativistic constituent qu- 
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ark model [6] can be viewed as an effective quantum field theory approach 
based on an interaction Lagrangian of hadrons interacting with their con- 
stituent quarks. Universal and reliable predictions for exclusive processes in- 
volving both mesons composed from a quark and antiquark and baryons com- 
posed from three quarks result from this approach. The coupling strength of 
hadrons H to their constituent quarks is determined by the compositeness con- 
dition Zh = [7,8] where Zh is the wave function renormalization constant 
of the hadron. The quantity Zjj is the matrix element between a physical 
particle state and the corresponding bare state. The compositeness condition 
Z H = enables us to represent a bound state by introducing a hadronic field 
interacting with its constituents so that the renormalization factor is equal to 
zero. This does not mean that we can solve the QCD bound state equations but 
we are able to show that the condition Zh = provides an effective and self- 
consistent way to describe the coupling of the particle to its constituents. One 
starts with an effective interaction Lagrangian written down in terms of quark 
and hadron variables. Then, by using Feynman rules, the S'-matrix elements 
describing hadron-hadron interactions are given in terms of a set of quark 
diagrams. In particular, the compositeness condition enables one to avoid the 
double counting of quark and hadron degrees of freedom. This approach is 
self-consistent and all calculations of physical observables are straightforward. 
There is a small set of model parameters: the values of the constituent quark 
masses and the scale parameters that define the size of the distribution of the 
constituent quarks inside a given hadron. The shapes of the vertex functions 
and the quark propagators can in principle be determined from an analysis of 
the Bethe-Salpeter (Fadde'ev) and Dyson-Schwinger equations, respectively, 
as done e.g. in [9,10]. In the present paper we, however, choose a more phe- 
nomenological approach where the vertex functions are modelled by Gaussian 
forms and the quark propagators are given by local representations. We have 
demonstrated in our papers [6,11] that the relativistic constituent model is con- 
sistent with the heavy quark symmetry in the limit of infinite quark masses. 
We mention that the authors of [12] have developed a relativistic quark model 
approach for meson transitions which shows some similarities to our approach. 
They also use an effective heavy meson Lagrangian to describe the couplings 
of mesons to quarks. They use, however, point-like meson-quark interactions. 
Loop momenta are explicitly cut off at around 1 GeV in their approach [12]. In 
our approach we use momentum dependent meson-quark interactions which 
provide for an effective cut-off of the loop integration. 

We have elaborated the so-called relativistic three-quark model (RTQM) to 
study the properties of heavy baryons containing a single heavy quark (bot- 
tom or charm). For the heavy quarks we have used propagators appropriate 
for the heavy quark limit. Various observables describing semileptonic and 
nonleptonic decays as well as one-pion and one-photon transitions have been 
successfully described in this approach [11]. Recently, the RTQM was extended 
to include the effects of finite quark masses [13]. 



In this paper we employ the RTQM [13] to calculate the form factors and 
widths of the semileptonic decays of the lowest lying S^ c and S cc baryons. We 
follow the strategy adopted in Ref. [5] where we have studied leptonic and 
semileptonic decays of the 5 c -meson. We employ the impulse approximation 
in calculating the matrix elements which previously has been widely used 
in phenomenological Dyson-Schwinger equation studies (see, e.g. Ref. [10]). In 
the impulse approximation one assumes that the vertex functions depend only 
on the loop momentum flowing through the vertex. 

We start with a brief description of our approach. As was mentioned above, 
baryons are described as bound states of constituent quarks in the RTQM. 
The general form of the SU(5)-invariant lagrangian describing the interaction 
of three low-lying SU(5)-multiplets with their three-quark currents are written 

as 

'MntV^J = *-int \ x ) + *-int \ x ) + Mnt \ X ) {*■) 

where 

C\£~ (x)=g F F [mim * ms] (x) J™ im2m *(x) + h.c, 
Cl^ + (x)=g B B [mim ^ ms (x) J™ im2m3 (x) + h.c, 
C-^ + (x)=g D D^ mim2m ^(x) J™ im2m ^(x) + h.c. 

Here m f = u, d, c, s, b are flavor indeces. According to the SU(5)-classification 

5 ® 5 © 5 = 10a © 40 M © 40 M © 35 5 

there is the antisymmetric decuplet F^ mim2m '^ with J p = | , two 40-plets 
^[m 2 m.i]m 3 w [^ m j X ed symmetry and J p = | and a symmetric 35-plet 
D { mim2m3 } with jP = |+_ 

The three-quark currents are written as 

Jp im2ni3 (x)= I dxi dx 2 dx3$ F (x;x 1 ,X2,x 3 ) 

x7 M 7 5 C(*i) (C(*2)cWcW 
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Table 1 

The lowest-lying states of double heavy H-baryons with diquarks in the symmetric 

state. The light quark q denotes u or d. 



Baryon 


Interpolating current 


Mass (GeV) 


^cs 


-V2 7 Vc a {s b C^q c )s abc 


2.47 


^cc 


7V<Za {c b C^c c )e abc 


3.61 


Z'bs 


-V2^ 5 b a (s b C^q c )e abc 


5.80 


^fec 


-^/2Yl b b a {c b C^q c )e ahc 


7.00 



where a^ = 1, 2, 3 are color indices. F, B and D denote the above three multi- 
plets. Note that the function $# is taken to be invariant under the translation 
x — > x + a which guarantees Lorentz invariance for the interaction Lagrangian 
Eq. (1). SU(5)-symmetry is broken by employing explicit baryon and quark 
mass values when calculating matrix elements. In this paper we limit our at- 
tention to the basic semileptonic decay modes of the lowest lying l/2 + double 
heavy baryons: S bcq ->• 2 bsq + Iv, H bcq ->• 2 ccq + Iv and H ccq ->• H csq + Iv (q=u 
or d). The appropriate interpolating currents with diquarks in the symmetric 
state and masses are shown in Table 1. The values of the masses are taken 
from potential models (see, for example, [2,14]). 

The semileptonic transition amplitude is defined as 



A(E i (p)^E f (p')lv 



v if ^{io»v) (*f(j/)K+f(p,rt»i(p) 



(2) 



where M = 7 M (1 — 7 5 )- V%f is the relevant element of the CKM- matrix where 
we use V bc = 0.04 and V cs = 0.97. The amplitude A M is decomposed into a set 
of six invariant form factors which are functions of the momentum transfer 
squared q 2 = (p — p') 2 only: 



7 M (pV _ pA 7 5) + ia p, q v ( F V _ pA 7 5) + g M [pV _ pA yj) 



(3) 



We shall not write down rate expressions in terms of these form factors since 
these have been worked out in great detail in Ref. [3]. 

In the impulse approximation which is being employed in our approach, the 
matrix element A M is calculated according to 
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A '-M""> = - 12 W J2^lj (££*<-* '>*/<-* '><*•/• 



(4) 



where 



OUc = -V2 7 a 7 5 S q (k 2 ) ^ S c (h + k 2 ) r S c (h + p') O" S b (h + p) 7V , 
CL+bs = 7°7 5 S b (h + p) 7V tr (S s {k 2 - q) 0» S c (k 2 ) ^ S q (h + k 2 ) i a ) , 
C^ cs = -v^ 7V S c (h + p) r S c (k 2 ) 0% S s (k 2 + qh 13 S q (k 2 - h) 7 a 7 5 , 

and where 0+ = 7 M (l+7 5 ) and k 2 = k 2 + {ki + k 2 ) 2 + k 2 . The quark propagator 
is chosen to have a local form 

s *(k) = u (i = u,d,s,c,b) (5) 

mi— f 

with mi being a constituent quark mass. The vertex function <f>H is directly 
related to the Fourier-transform of the function $# 



' / j &iPi 



$ H (pi,...,p 4 )=Jdx 1 ...Jdx 4 e *=i $ H (x u ...,x 4 ) 



i=i 



Generally, <pH is a function of three momentum variables. However, in the 
impulse approximation employed in our approach, we assume that it only 
depends on the sum of relative momentum squared as indicated in Eq. (4). 

The compositeness condition reads 

Z H = l- g 2 H Y>'{m H ) = (6) 

where E'(mn) is the derivative of baryon mass operator taken on its mass- 
shell. In the impulse approximation Eq. (6) may be rewritten in a form suitable 
for the determination of the coupling constants: 



2 / a m / a fc 2 , 2 ,_, 2 \ ™ I _ ix (7) 

^VqWiqz / (c>„\\; / (r>~\4,; ^919293 V ^ J -^919293 I ^"fl ' ' VV 



rf 4 /ci /■ rf 4 /c 2 
(2tt) 4 U (2vr) 4 i 

^9293 = 7°7 5 5 91 (fc + p) 7 M 5 W (h + p) 7 V tr (S g2 (h + k 2 ) 1 a S q3 (k 2 ) /) , 
(919293) = ( fo cg) , (6sg) , (csg) , 

£>£, = 7 a 7 5 S 9 (* 2 ) 7 /3 7 5 tr ( 7 a ^(fei + p) 7" S c (h + p) ^ S c (h + k 2 )) . 
For the coupling constants we obtain g bcq = 0.96, g bsq = 3.33, g ccq = 2.63, 

gcsq = <J- '"■ 



Next we turn to the calculation of the transition form factor. The calculational 
techniques are outlined in Ref. [13]. The three main ingredients are 



use of the Laplace transform of the vertex function 

oo 

$(z) = fds<$> L (s)e~ sz 



the a-transform of the denominator 

1 



m? — (k + p) 2 



dae~ a(m ' 2 ' {k+p)2) 



• differential representation of the numerator 

(m+ J6+ f) e kq = I ™ + ^—+ A e kq 

The calculation of the transition form factors amounts to a two-loop integra- 
tions. Four of the eight two-loop integrations are done analytically. One ends 
up with 4-fold integrals which are not difficult to evaluate numerically. All 
calculations are done by using computer programs written in FORM for the 
manipulations of Dirac matrices and in FORTRAN for numerical evaluations. 

The common structure of the expressions for the form factors may be written 

as 



v(q 2 



J dtt3 Jw 6 i 1 ~ E *) { F{Z)W ° ~ \ F ^ Wl + \ F *W W *) ( 8 ) 



where 



F(z) = <f>in(z)<f>out(z) , Fi(z) = / drr i F{z + r) , 



o 



t Yl a i m \ ~ a iP 2 ~ a 2P' 2 + t 2 A 1 l(a 1 p - a 2 p 



z „ . /L ^^. l .,. q 
and where the matrix A is defined as 



/\2 



2 + t(ai + a 2 + a 3 ) 1 + to 3 

1 + ta 3 2 + t(a 3 + a 4 ) 



The functions Wi contain integration variables and masses. One has to empha- 
size that the above expressions are valid for any vertex functions decreasing 



rapidly enough in the Euclidean region. Since the quark masses satisfy the 

3 

confinement constraint tuh < J2 m qi all form factors are real. 

Before presenting our numerical results we need to specify our values for the 
constituent quark masses and shapes of the vertex functions. As concerns the 
vertex functions, we found a good description of various physical quantities 
[5,6,11,13] adopting a Gaussian form. Here we apply the same procedure using 
4>H(k 2 E ) = exp{— k 2 E / A 2 H } in the Euclidean region. The magnitude of A H char- 
acterizes the size of the vertex function and is an adjustable parameter in our 
model. The A# parameters in the meson sector were determined [6] by a least- 
squares fit to experimental data and lattice determinations. The nucleon Ajv 
parameter was determined from the best description of the electromagnetic 
properties of the nucleon [6] . The Ah parameters for baryons with one heavy 
quark (bottom or charm) were determined by analyzing available experimen- 
tal data on bottom and charm baryon decays. Since there is no experimental 
information on the properties of double heavy baryons yet we use the sim- 
ple observation that the magnitude of Ah is increasing with the mass value 
of the hadron whose shape it determines. Keeping in mind that Ajy = 1-25 
GeV, Aq qq = 1.8 GeV and A Bc = 2.43 GeV, we simply choose the value of 
A<3Q g = 2.5 GeV for the time being. We found that variations of this value by 
10 % does not much affect the values of form factors. We employ the same 
values for the quark masses (see, Eq.(9)) as have been used previously for the 
description of light and heavy baryons [6,11]. We thus use 



m u m s m c m b 

(9) 



0.420 0.570 1.67 5.06 (GeV) 



The resulting form factors are approximated by the interpolating form 

f(i 2 ) = i — f -Pr—, (io) 

1 - a x q l + a 2 g 4 



It is interesting that for most of the form factors the numerical fit values of 
a\ and a<i obtained from the interpolating form (10) are such that the form 
factors can be represented by dipole formula 



J(0) 
q 2 /m 2 



f{q 2 ) * d{q 2 ) = j, J V> . , (11) 



The values of my in the dipole representation are very close to the values 
of the appropriate lower- lying (qq 1 ) vector mesons (m£>.=2.11 GeV for (c-s)- 
transitions and mg* ~ mg c =6.4 GeV for (b-c)-transitions). In Fig.l we show 



two representative form factors and their dipole approximations. We have 
shown in Ref. [5] that the form factors of the CKM-enhanced semileptonic 
.B c -decays may be approximated by a monopole function. It is gratifying to 
see that our relativistic quark model with the Gaussian vertex function and 
free quark propagators reproduces the monopole in the meson case and the 
dipole in the baryon case for most of the form factors. 

Finally, in Table 2 we present our predictions for the decay rates and compare 
them with the free quark decay width which is the leading contribution to the 
semileptonic inclusive width (x = itiq Jfn 2 Q x ) 



r (^/) = |%| 2 i|^f (l-8x + 8x 3 -x 4 -12x 2 lnx) . (12) 



The nonleading corrections to the leading order rate Eq.(12) lower the inclusive 
rate by approximately 1% and 15% in the b — ► c and c — ► s case, respectively 
(see e.g. [15]). In the numerical evaluation of the inclusive rate Eq.(12) we used 
current pole mass values m& = 4.8 GeV and m c = 1.325 GeV [16]. There are 
many values quoted in the literature for the current pole mass of the strange 
quark. For the sake of definiteness we take m s = 0.15 GeV. Note that in both 
cc — > cs and bc — ► cc decays there is an additional factor of 2 due to the 
fact that there are two c-quarks in the double charmed baryon in the initial 
and final state, respectively [2,4]. One notes that the rates for the exclusive 
modes H& c — > S cc + Iv and H& c — » H& s + Iv are rather small when compared 
to the total semileptonic inclusive rate estimated by Eq.(12). The remaining 
part of the inclusive rate would have to be filled in by decays into excited or 
multi-body baryonic states. Note that the smallness of the exclusive/inclusive 
ratio of the above exclusive modes markedly differs from that of the mesonic 
semileptonic b — > c transitions, where the exclusive transitions to the ground 
state S-wave mesons B — > D, D* make up approximately 66% of the total 
semileptonic B — > X c rate [17]. For A& — > A c transitions one expects even 
higher semileptonic exclusive-inclusive ratio of amount 80% [15]. Note that 
the rate for Sfe c — > S cc + Iv is of the same order of magnitude as the rates 
calculated for the corresponding double heavy mesonic decays B c — > f] c + Iv 
and B c — * J/ty + Iv [5]. The QCD sum rule and potential model predictions 
for the rates of St c — ► S cc + Iv and St c — ► E cs + Iv given in [2] exceed our 
rate predictions by factors of 10 and 3 respectively. In fact, the exclusive 
semileptonic rates given in [2] tend to saturate the inclusive semileptonic rates 
calculated from Eq.(12) as given in Table 2. 

In Table 3 we present values for the invariant form factors at g^ in = and 
9max = {rrii — rrif) 2 . Note that the values of the axial vector form factor Ff 
are rather small for the two decays E& c — > S cc + Iv and S^ c — > Sb s + Iv. This 
provides for a partial explanation of why the rates of these two modes are small 
compared to the inclusive semileptonic rate. Also the zero recoil values of the 



vector form factors are significantly below the value of one which one would 
expect from a naive application of the heavy quark limit. The smallness of 
the vector form factors provide for the remaining explanation of the smallness 
of the predicted respective rates. We mention that the QCD sum rule and 
potential model estimates of the zero recoil values of both the vector and 
axial vector form factors FY and F^ given in [2] are close to one. In the 
model of [2] the form factors _F 2 y and F^ are set to zero. In our approach 
we find that the numerical values of F^ and F^ are quite small compared to 
those of FY and F^ in all cases when expressed in terms of the mass scale 
(rrii + rrif). 

It is interesting to compare our full zero recoil results with those of a naive 
spectator quark model calculation where the zero recoil values of the form 
factors result from a simple overlap calculation. For the spin-flavor wave func- 
tions we use the SU (12) spin-flavor wave functions of the naive spectator 
quark model (6 flavors x 2 spin projections) as extended from its original 
SU(6) version to include the heavy flavors. As mentioned before this does not 
imply that one is assuming SU (12) symmetry to hold for the transitions since 
one is using physical quark and baryon masses which badly break the SU (12) 
symmetry. Instead one uses SU (12) symmetry only to construct the spin-flavor 
wave functions of the heavy and double heavy baryons. In the SU (12) naive 
spectator quark model the spin-flavor wave functions \B > of S& cg , S{, sg and 
S ccg baryons with diquarks in the symmetric state and positive (+1/2) spin 
projection are given by 



1 

6^2' 
1 

6^2' 
1 

2^2 



z bcq ; t> = — y=\2qcb + 2cqb - qbc - cbq - bqc - bcq > \ TIT + ITT "2 TTI> 
l bsq ; t> -= iri =\2qsb + 2sqb - qbs - sbq - bqs - bsq > | TIT + ITT -2 TTI> 
i ccq ; |> = ^^Igcc + cqc - 2ccq > | TIT + 1 11 "2 TTI> 



The values of FY and Ff form factors at zero recoil can be calculated from 
the matrix elements 



FY =< B'\ Y.\Ifi\ {i) \B > and F, A =< B'\ J2<T 3 I fl ]®\B > 

i=\ i=l 

where 03 is the Pauli (third component) spin matrix and Ifi is the flavor 
matrix responsible for the s.l. transitions. 

In the Table 4 we present the results for the FY and Ff form factors at 
zero recoil calculated in the naive spectator quark model. It is evident that 
the spectator quark model prediction of a vanishing axial vector coupling 



F 1 A (0) = for the decay E bc — > E cc + lu is close to the suppressed value of 
F{ 4 (0) = —0.091 in the full calculation. The ratios of the axial and vector 
couplings in the decay E bc — > E bs + lu are F 1 A (0)/F 1 y (0) = 1/2 (naive quark 
model) and F^(0)/F^(0) ps 1/6 (our approach). For S cc — »■ S cs + /z/ decay one 
has F{ 4 (0)/F 1 y (0) = 1 (naive spectator quark model) and F{ 4 (0)/F 1 y (0) « 1.3 
(our approach). The different normalizations of axial and vector constants 
obtained in our approach compared to the naive spectator quark model result 
and the suppression of the ratio F 1 A (0)/F 1 v '(0) in the decay E bc — > E bs + lu 
in the full calculation can be explained by relativistic effects and nontrivial 
heavy quark/baryon mass dependence of the relevant matrix elements. 

In order to test the sensitivity of our results on various choices of the bottom 
and charm quark masses, we have varied their values within a reasonable 
range. From the confinement constraints one obtains lower permissible values 
for heavy quark masses: m c > (m= cc — m u )/2 = 1.60 GeV and m b > (m Sbc — 
m c — m u ) = 4.98 GeV. The upper values were found from the experimental 
bounds for the A& — > A c e~u and A c — > A s e + u decay rates: m c < 1.72 GeV 
and m b < 5.25 GeV [13]. We have calculated the values of decay rates of 
double baryons for three set of quark masses in Table 5. The decay rates do 
not change significantly in the chosen regions of the heavy quark masses. 
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Russian Fund of Basic Research 01-02-17200 and the DFG (FA67/25-1). 
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Table 2 

Calculated decay widths of lowest lying J p = l/2 + double heavy H-baryons. Inclu- 
sive widths are calculated using the current quark pole masses. 



Mode 


Decay widths, ps * 


RTQM 


Inclusive width 


u be ► il cc + I V 


0.012 


2-r (6^c) = 0.162 


2 be — *• S bs + ^ 


0.043 


r (c^s)= 0.122 


Ci cc ► s cs + ly 


0.224 


2-r (c^s)=0.244 



Table 3 

Values of F^ and F^ form factors at maximum and zero recoil. 





Z-bc ~ 


"* ^cc 


-6c " 


-*■ Sft s 


^cc 


■* 1 — 'CS 


FY 


Ff 


FY 


^i 4 


FY 


Ff* 


<f 


= 


0.46 


-0.091 


0.39 


0.061 


0.47 


0.61 


<f 


= Q 2 


0.83 


-0.086 


0.58 


0.065 


0.59 


0.77 



Table 4 

Values of FY and FT form factors at zero recoil in the naive spectator quark model. 



Quantity 


^beq > ^ccq 


^beq * ■— 'bsq 


^ccq * ^csq 


FY 


1/V2 


1 


1/V2 


Ff 





1/2 


1/V2 



Table 5 

Calculated decay widths of double heavy H-baryons for three sets of quark masses. 



Decay 


m c = 1.60 GeV 
m b = 4.98 GeV 


m c = 1.67 GeV 
m b = 5.06 GeV 


m c = 1.72 GeV 
m b = 5.25 GeV 


T(H bc -»■ ~ CC + /P), ps -1 


0.0102 


0.0117 


0.0123 


r(H bc -»■ E hs + lP), ps -1 


0.0498 


0.0432 


0.0371 


r(H cc ->■ E cs + lv), ps -1 


0.258 


0.224 


0.208 
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Fig. 1. Upper panel: Form factor F± (q 2 ) (solid-dotted line) for be — ► bs transitions 
and its dipole approximation (solid line) with m cs = 2.88 GeV; Lower panel: Form 
factor FY (q 2 ) (solid-dotted line) for be — > cc transitions and its dipole approximation 
(solid line) with m^ = 6.81 GeV. 
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